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Abstract 

We consider a wide class of semi linear Hamiltonian partial differential equa- 
tions and their approximation by time splitting methods. We assume that the 
nonlinearity is polynomial, and that the numerical trajectory remains at least uni- 
formly integrable with respect to an eigenbasis of the linear operator (typically 
the Fourier basis). We show the existence of a modified interpolated Hamiltonian 
equation whose exact solution coincides with the discrete flow at each time step 
over a long time depending on a non resonance condition satisfled by the stepsize. 
We introduce a class of modified splitting schemes fulfilling this condition at a 
high order and prove for them that the numerical fiow and the continuous fiow 
remain close over exponentially long time with respect to the step size. For stan- 
dard splitting or implicit-explicit scheme, such a backward error analysis result 
holds true on a time depending on a cut-off condition in the high frequencies 
(CFL condition). This analysis is valid in the case where the linear operator has 
a discrete (bounded domain) or continuous (the whole space) spectrum. 

MSG numbers: 65P10, 37M15 

Keywords: Hamiltonian interpolation, Backward error analysis. Splitting inte- 
grators, Non linear Schrodinger equation, Non linear wave equation. Long-time 
behavior. 



1 



1 Introduction 



The Hamiltonian interpolation of a symplectic map which is a perturbation of 
the identity is a central problem both in the study of Hamiltonian systems and 
in their discretization by numerical methods. This question actually goes back to 
Moser [23] who interpreted such a map as the exact flow of a Hamiltonian system 
in the finite dimensional context. Such a result was later refined and extended by 
Benettin and Giorgilli [3] to the analysis of symplectic numerical methods, and 
leads to the seminal backward error analysis results of Hairer, Lubich [16] and 
Reich [24] for numerical integrators applied to ordinary differential equations. 
These results constitute now a cornerstone of the modern geometric numerical 
integration theory (see [18, 25]). 

In the finite dimensional case, the situation can be described as follows: if 
{p,q) I— > ^^{p,q) is a symplectic map from the phase space M^*^ into itself, and 
if this map is a perturbation of the identity, ~ Id + 0{h), then there exists 
a Hamiltonian function Hfi{p,q) such that "i>^ can be interpreted as the flow at 
time t = h of the Hamiltonian system associated with Hh- If is analytic, 
such results holds up to an error that is exponentially small with respect to 
the small parameter h, and for {p, q) in a compact set of the phase space. In 
application to numerical analysis, the map ^^{p,q) ~ '^^{p^q) is a numerical 
approximation of the exact flow associated with a given Hamiltonian H. As a 
consequence, the modified Hamiltonian Hh-, which turns out to be a perturbation 
of the initial Hamiltonian if, is preserved along the numerical solution, q-n) = 
(^''^)"'(po, 9o), n G N. More precisely we have 

{Pn.qn) = {■^^T{po,qo) = m\)iPo,qo)+nheM-'^/ich)) (1.1) 

under the a priori assumption that the sequence {pn, '?n)neN remains in a compact 
set used to derive the analytic estimates. Here the constant c depends on the 
eigenvalues of the quadratic part of the Hamiltonian H (i.e. the linear part of 
the associated ODE). As a consequence the qualitative behavior of the discrete 
dynamics associated with the map over exponentially long time can be pretty 
well understood through the analysis of the continuous system associated with 
Hh. 

The extension of such results to Hamiltonian partial differential equations 
(PDE) faces the principal difficulty that the Hamiltonian function involves oper- 
ators with unbounded eigenvalues making the constant c in the previous estimate 
blow up. The goal of this work is to overcome this difficulty and to give Hamilto- 
nian interpolation results for splitting methods applied to semi linear Hamiltonian 
PDEs. 

We consider a class of Hamiltonian PDEs associated with a Hamiltonian H 
than can be split into a quadratic functional associated with an unbounded linear 
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operator, i^o, and a nonlinearity P which is a polynomial functional at least cubic: 

H = Ho + P. 

Typical examples are given by the non linear Schrodinger equation (NLS) 

idtu = Au + f{u,u) (1.2) 

or the non linear wave equation (NLW) 

dttu-Au = g{u) (1.3) 

both set on the torus T'^ or on the whole space M'^. Here / and g are polynomials 
having a zero of order at least three at the origin, for instance f{u,u) = \u\'^u for 
the cubic defocusing (NLS) and g{u) = v? for the classical (NLW). 
To approximate such equations, splitting methods are widely used: They consist 
in decomposing the exact flow at a small time step h as compositions of the 
flows and <I>p. The Lie- Trotter splitting methods 

u{nh,x)^Un{x):={^'ho°^pTM^) or n„(x) := ($^ o $^J"no(x) (1.4) 

are known to be order 1 approximation scheme in time, when the solution is 
smooth. Here smooth means that the numerical solution belongs to some Sobolev 
space , s > 1, uniformly in time, see [21, 19] for the linear case and [22] for 
the non linear Schrodinger equation. 

These schemes are all symplectic, preserve the norm if Hq and P do, and 
can be easily implemented in practice. For instance, in the cases of (NLS) or 
(NLW) just above, we can diagonalize the linear part and integrate it by using 
fast Fourier transform and integrate the nonlinear part explicitly (it is an ordinary 
differential equation). More generally, instead of fast Fourier transform, pseudo 
spectral methods are used if the spectrum of is known and available. 

The Hamiltonian interpolation problem described above can be formulated 
here as follows: is it possible to find a modified energy (or modified Hamiltonian 
function) Hh depending on Hq, P and on the chosen stepsize, such that 

^p°'&Ho=^'J>k- (1-5) 

Formally, this question corresponds to the classical Baker-Campbell-Hausdorff 
(BCH) formula (see [2, 20]) for which the modified Hamiltonian Hh expresses as 
iterated Poisson brackets between the two Hamiltonian P and Hq. Hence we ob- 
serve that the validity of such representation a priori depends on the smoothness 
of the discrete solution Un- But this is not fair, as there is no reason for u„ to be, 
a priori, uniformly smooth over long time. Actually, it is well known that this 
assumption is not satisfied in general - even if the exact solution of the continuous 
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system is smooth - making such analysis irrelevant in many cases of applications, 
see for instance [10]. Note that this smoothness assumption balances the blow up 
of the constant c in the finite dimensional exponential estimates (1.1). In other 
words, the norm for which (1.5) is valid with the "finite dimensional" method 
is not the norm required for the analytic estimates of the error in the infinite 
dimensional context. 

In this work, we give backward error analysis results in the spirit of [3, 16, 24] 
for the map (1.5). We follow a general approach recently developed in [8] for 
linear PDEs. In this linear context, Debussche and Faou proved a formula of 
the form (1.5) by considering smoothed schemes, namely schemes where ^^^^ is 
replaced by the midpoint approximation of the unbounded part Hq. This yields 
schemes of the form 

<^'ho^\^ or 1>ii„o$^ (1.6) 
where Aq depends on h, and is an approximation of the operator hHo when 

In the present work we consider the context of nonlinear PDEs [P is at least 
cubic) and we assume that the eigenvalues of Aq are of the form 

Aa = OLhihiOa) (1.7) 

where we denote by oJa the frequencies of Hq (here the index a belongs to a discrete 
set or is a continuous variable), and where is a ^/ier function satisfying ah{x) ~ 
X for small x. It turns out that under such an assumption, (1.6) remains an 
order one approximation of the continuous solution, provided that the numerical 
solution is smooth. This general setting contains the standard splitting methods 
(1.4) which corresponds to the choice ah{x) = x for all x, and the midpoint 
approximation of the linear flow which corresponds to the function (see [8]) 

ah{x) = 2 arctan(x/2). 

Our results can be described as follows: we assume that, in the Fourier vari- 
ables associated with the diagonalization in an eigenbasis of Hq, the numerical 
trajectory remains bounded in the space i]., for some fixed s > 0. In the case of 
the Laplace operator on T*^, this means to assume that Un has its Fourier trans- 
form satisfying {\kYu^{k)j.^id) in (^{lA) while in the case of the same operator 
on M'^, this corresponds to assume that u„ has its Fourier transform satisfying 
|^|*?x^(^) in L^(R'^). When for instance s = 0, this implies in both cases that the 
numerical trajectory remains bounded in the L°° space associated with the space 
variable uniformly in time. This is actually the standard assumption is the finite 
dimensional case: no blow up in finite time. 



4 



Then for z (the Fourier variable) in a fixed ball of , we construct a Hamiltonian 
Hh such that 

\\^\ o ^\^z) - < Wll,. < /i^+^(CiV)^ (1.8) 

for some constant C depending on P and s. It turns out that the number 
depends on a sort of non resonance condition satisfied by the eigenvalues of the 
operator Aq. With Aq the eigenvalues of Aq (see (1.7)) this condition can be 
written 

Vj = l,...,r, V(ai,...,aj) |Aai±---±AaJ < 27r. (1.9) 
If this is satisfied for some r > 2, then we can take in the previous estimate 

N ={r- 2)/{rQ - 2) 

where tq is the degree of the polynomial P. Note that this condition is indepen- 
dent of the regularity index s. We thus see that the actual existence of a modified 
PDEs almost coinciding with the numerical scheme at each time step depends on 
the degree of the polynomial P and on the non resonance condition (1.9). 
For numerical schemes associated with a filter function 

ah{x) = arctan(x/\//i), (1-10) 

we can actually prove that the non resonance condition (1.9) is satisfied for large 
r ~ Xjyfh. The analytic estimate (1.8) then yields an exponentially small error at 
each step. The choice (1.10) typically induces a stronger regularization in the high 
frequencies than the midpoint rule, without breaking the order of approximation 
of the method. 

For the other classical schemes (standard splitting, implicit-explicit schemes), 
the non resonance condition (1.9) is in general not satisfied unless a Courant- 
Friedrichs-Lewy (CFL) condition is imposed [7], depending on the wished number 
r. Although the CFL number tends to when r goes to oo for all schemes, we 
obtain very different results for relatively small r regarding on the method and on 
the specific form of the nonlinearity. For example the combination of the midpoint 
approximation of the linear part in an implicit-explicit scheme together with a 
CFL condition with CFL number of order 1 in the case of the cubic nonlinear 
Schrodinger equation yields a modified equation with precision of order hJ at 
each step (see section 5 for details and other examples). 

We can precise in which sense the estimate (1.8) induces a better control of 
the numerical solution. Actually, the modified Hamiltonian depends on and 
P and can roughly speaking be described as a Hamiltonian A^^jh^P where P is 
a modified nonlinear Hamiltonian. For example in the case of the method (1.10) 
applied to the cubic (NLS) on M"^ we show (cf. section 6) that the modified energy 
reads 

/ ^arctan(^/^|eP)|n(Ol'de + ^'i'^M, (I-H) 
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and is preserved up to an error 0{h) over exponentially long times t < e^^^, 
provided the numerical solution remains bounded in £^ = only. Here, u denotes 
the Fourier transform of u and pj^^\u) is a polynomial nonlinearity of degree 4. 

The actual expression of P^^^(u) depends on small denominators controlled by 
the non resonance condition (4.1) and is given in Theorem 6.5. Note that this 
modified energy will be close to the original energy only if u is smooth, something 
that is not guaranteed by the present analysis unless the a priori regularity index 
s is large enough. However, the preservation of (1-11) implies the control of the 
norm of low modes of the numerical solution, and of the norm of high 
modes, as in [8], over very long time. 

The relative question of persistence of smoothness of the numerical solution 
has recently known many progresses: see [9, 11, 12, 17, 6, 13, 14] and [10]. How- 
ever, we emphasize that the resonances considered in all these works are not of 
the same type as here. They read 

Aai ± • • • ± Aa^ ~ 2A;7r for some /c G Z (1.12) 

and they control the persistence of the regularity of the solution over long time. 
As explain in [11, 12, 10], the case k = corresponds to physical resonances (they 
actually are close to resonances of the PDE under study, see also [4, 15]) while 
the case k ^ corresponds to numerical resonances (they are not resonances of 
the PDE). 

In the present work, we address a different question and we point out that possible 
physical resonances do not contradict the existence of the modified energy (1.11). 
Such relations are in fact also resonances of the modified equation associated with 
the modified Hamiltonian Hh . In such a case, no preservation of the regularity in 
high index Sobolev spaces H"^ for the numerical solution can be expected, though 
the existence of the modified equation is guaranteed by the present analysis even 
if s = 0. Only numerical resonances can avoid the existence of a modified energy. 

We end this introduction by two important commentaries: 

Commentary 1.1 We emphasize that we use the assumption that the nonlin- 
earity has zero momentum. This ensures the fact that the Poisson bracket of 
two Hamiltonian with bounded polynomial coefficients will still act on the space 
i], and that the corresponding vector field will be well defined on i]. This kind 
of assumption is made in [12] and precludes Hamiltonian PDEs with strong non 
local nonlinearity. Note however that we may obtain results with small momen- 
tum by combining the results in [8], i.e. by assuming that the nonlinearity P{x, u) 
depends smoothly on x and polynomially on u. 

Commentary 1.2 Note that similar results can be obtained for Strang-like split- 
ting methods generalizing (1.6), i.e. schemes of the form 
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for which the analysis is very close to the one developed here. Similarly, we 
conjecture that the same kind of results holds true for any symplectic Runge- 
Kutta method under CFL condition, but this leads to technical difficulties. The 
big advantage of the splitting approach is that the induction equations building 
the modified Hamiltonian directly express at the level of the Hamiltonian, and not 
on the corresponding vector field (in particular we do not have to show that the 
vector field constructed at each step is - indeed - Hamiltonian). 

2 Abstract Hamiltonian formalism 

In this section we describe the general form of Hamiltonian PDEs we are consid- 
ering. We focus on the case where the spectrum of the linear operator is discrete: 
elliptic operator on the torus or on a bounded domain with boundary condi- 
tions. We give in the last section of this paper the principal changes that have 
to be made to obtain similar results in the case where the spectrum of the linear 
operator is continuous (Schrodinger or wave equation over the whole space). 

So typically, the solution of the PDE under study is decomposed in the eigen- 
basis of the linear part: 

ii^{t,x) = Y,ik{t)(i>k{x) 

and we observe the PDE in the Fourier-like variables ^ = {^k)- 

The main difference with the presentation in [11, 12, 15] lies in the choice 
of the phase space. We consider Fourier variables ^ belonging to l\ and not 
to ll the weighted I space which corresponds to function -0 belonging to the 
Sobolev spaces . In particular, considering Fourier coefficients in £j remains 
to consider bounded functions, i.e. to assume that the solution (or its numerical 
approximation) of the PDE under study is (essentially) bounded. 

2.1 Setting and notations 

We set M = or N*^ (depending on the concrete application) for some d > 1. 
For a = (oi, . . . , as) G N , we set 

\a\ = max (1, -|- • • • -|- a^j . 

We consider the set of variables (^a, rjh) £ C-'^ x C-^ equipped with the symplectic 
structure 

i^d^a/\drja. (2.1) 

We define the set Z = M x {±1}. For j = (a, 6) G Z, we define \ j\ = \a\ and we 
denote by j the index (a, —6). 
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We will identify a couple r?) G C-^ x C-^ with {zj)j(zz G via the formula 
By a slight abuse of notation, we often write z = {^,7]) to denote such an element. 



r]a if S = -1. 



Example 2.1 In the case where Hq = —A on the torus T'^, the eigenbasis is the 
Fourier basis, N ='L'^ and ^ is the sequence associated with a function ip while rj 
is the Fourier sequence associated with a function (p via the formula 

For a given s > 0, we consider the Banach space il := C) made of elements 
z G such that 

Ikll^i := X] l-^'I'l^jl < 

and equipped with the symplectic form (2.1). We will often write simply £^ = ^o- 
We moreover define for s > 1 the Sobolev norms 



1/2 

For a function F of C^{il, C), we define its gradient by 



where by definition, we set for j = {a,6) £ M x {±1}, 

dF 

dF_ 

dzj 

Let F[{z) be a function defined on il. If H is smooth enough, we can associate 
with this function the Hamiltonian vector field Xniz) defined by 

Xh{z) = JVH{z) 

where J is the symplectic operator on il induced by the symplectic form (2.1). 
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For two functions F and G, the Poisson Bracket is (formally) defined as 

We say that z £ is real when = for any j £ Z. In this case, we write 
z = (C,^) for some ^ G C-^. Further we say that a Hamiltonian function H is 
real if /^(-z) is real for all real z. 

Example 2.2 Following Example 2.1, a real z corresponds to the relation (p = 
and a typical real Hamiltonian reads H{z) = Jj^ h{tpix), 4>ix))dx where h is a 
regular function from to C satisfying h{(, C) G M for all C, £ C 

Definition 2.3 For a given s > 0, we denote by Tig the space of real Hamiltoni- 
ans P satisfying 

P£C^{£l,C), and Xp £ C^{il,£l). 
Notice that for F and G in 7is the formula (2.2) is well defined. 

Remark 2.4 At this stage we have to note that, in general, the quadratic Hamil- 
tonian Hq corresponding to the linear part of our PDF will not be defined for 
£ £l. For instance for the NLS equation, Hq is associated with the Laplace 
operator (see example 2.1 above) and reads in Fourier variable ^'^\^k\'^ ■ Never- 
theless it generates a flow which maps £] into £] and which is given for all time 
t and for all indices k by ^fc(t) = e~*^*^fc(0). 

We will verify later (see example 2.7) that the typical real Hamiltonians given 
in Example 2.2 belong to the class Tis- Actually the proof is not totally trivial 
because the Fourier transform is not well adapted to the £] space. 

With a given Hamiltonian function H £ 7is, we associate the Hamiltonian 
system 

z = JVH{z) 

which can be written 

dH 

dH (2-3) 

Va = iT^(?,^) a£Af. 
OU 

In this situation, we define the flow ^^fji^) associated with the previous system 
(for an interval of times t > depending a priori on the initial condition z). Note 
that if z = (C,^) and if H is real, the flow (C*j^*) = ^*h^^) ^^^^ '^^^l fo'^ 
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time t where the flow is defined: ^* = r/*. When H is real, it may be useful to 
introduce the real variables pa and qa given by 

= -^iP<^ + and = -^(p^ _ iq^-^^ 

the system (2.3) is then equivalent to the system 

dqa 
dpa 

where by a slight abuse of notation we still denote the Hamiltonian with the same 
letter: H{q,p) = H{C,C). 

We now describe the hypothesis needed on the Hamiltonian nonlinearity P. 

Let £ > 2. We consider j = {ji, ■ ■ ■ , j^) G Z^, and we set for alH = 1, . . . ^ 
ji = (oj, 6i) where Ui £ M and 5i G {±1}- We define 

j = Ui,---,ji) with j^ = {ai,-6i), i = l,...J. 

We also use the notation 

^3 ~ ' ' ' '^k ■ 
We define the momentum A^(j) of the multi-index jf by 

M{i) = axbi^---^aibi. (2.4) 

We then define the set of indices with zero momentum 

Tt = {3 = {h.---.k)^Z\ with 7W(j) = 0}. (2.5) 

We can now define precisely the type of polynomial nonlinearities we consider: 

Definition 2.5 We say that a polynomial Hamiltonian P £ Vk if P is real, of 
degree k, have a zero of order at least 2 in z = 0, and if 

• P contains only monomials UjZj having zero momentum, i.e. such that 
-^(j) = when aj ^ and thus P formally reads 

k 

^(^) = E E (2-6) 
e=2 j&Ti 

with the relation aj = aj. 
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The coefficients aj are bounded, i.e. satisfy 

y£ = 2,...,k, Vj = (ji, • • • ,j» G X^, \aj\<C. 
In the following, we set 



p\\ = sup \aj\. (2.7) 



e=2 J'^^* 

Definition 2.6 We say that P G SVk if P £ Vk has coefficients aj such that 
Qj 7^ implies that j contains the same numbers of positive and negative indices: 

tJ{^ I ji = {ai, +1)} = ^{i I ji = {ai, -1)}. 

In other words, P contains only monomials with the same numbers of and rji. 
Note that this implies that k is even. 

Example 2.7 Following example 2.2, P{z) = JjdP{ip{x),(p{x))dx, where p is a 
polynomial of degree k in C[X,Y] satisfying p{C,C) £ ^ and having a zero of 
order at least 2 at the origin, defines a Hamiltonian in Vk ■ 
An example of polynomial Hamiltonian in SVik is given by P = f \tlj\'^^dx. 

The zero momentum assumption in Definition 2.5 is crucial in order to obtain 
the following Proposition: 

Proposition 2.8 Let k >2 and s > 0, then we have Vk C Us, and for P E Vk, 

we have the estimates 

P{z)\ < \\P\\ (^maxjlzli;) (2.8) 



■ n=2,...,fc 

and 



Vze4, \\Xp{z)\\^,<2k{k-iy\\P\\\\z\\^,{ me^i \\z\\l). (2.9) 
Moreover, for z and y in £l, we have 



\\Xpiz)-Xp{y)\\^, < Ak{k-ir\\P\\ { max (||?/||; , ||z||;)) . (2.10) 

S Tt J_ , . . . J A, s s s 

Eventually, for P G Vk and Q G Vi, then {P, Q} £ Vk+e~2 and we have the 
estimate 

\\{P,Q}\\ <2ke\\P\\ WQW . (2.11) 
// now P e SVk and Q £ SVk, then {P, Q} £ SVk+i-2- 
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Remark 2.9 The estimate (2.9) is a sort of tame estimate. Notice that the same 
estimate with £l replaced by (-1 is proved in [15] under a decreasing assumption on 
the coefficients of the polynomial P. Actually the present proof is much simpler. 
We also notice that, with this estimate and following [15], we could develop a 
Birkhoff normal form theory in £\ . The only disagreement with this choice of 
Fourier space is that is not the image by the Fourier transform of L°° , making 
difficult the pull back of the normal form to the original variables. 

Proof. Assume that P is given by (2.6), and denote by Pi the homogeneous 
component of degree i of P, i.e. 

Pi{z) = OjZj, i = 2,...,k. 

jeii 

We have for all z 

\Piiz)\ < \\Pi\\ \\z\t < \\Pi\\ \\z\\' 



1^1 



The first inequality (2.8) is then a consequence of the fact that 



\P\\ 



(2.12) 



i=2 



Now let j = (a, e) £ Z he fixed. The derivative of a given monomial Zj 
with respect to zj vanishes except if j C j. Assume for instance that j = ji. Then 
the zero momentum condition implies that . . . , ji^i) = —ea and we can 

write 

dP^ 



\J\ 



dzj 



<m\\ E (2-13) 

jeZi-\M{j)=-ea 

Now in this formula, for a fixed multiindex jf, the zero momentum condition 
implies that 

br < (Uil + • • • + < - 1)' max 



n=l,...,i— 1 



(2.14) 



Therefore, after summing in a and e we get 

\\Xp^{z)\\^,<2i{i-ir\\P,\\ E max | 



< 2i{i - iy\\Pi\ 



(2.15) 



which yields (2.9) after summing in z = 2, . . . , /c. 

Now for z and y in i], we have with the previous notations 

^ I ., o dPj 



dZjdZq 



[ty + (1 - t)z)dt 



yq\ 
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But we have for fixed j = (e, a) and q = {5, b) in and for all u £ il 
OR 



dZjdZq 



(u) 



b'l ■ ■ ■ %i-2i- 

jG2»-2 ^M{j)=-<ia-Sb 

In the previous sum, we necessarily have that M.{j,j,k) = 0, and hence 

i-2 

br < (ijii + • • • + \j^-2\ + \q\y <{^- ir kr n i-^^i'- 

n=l 

Let u(t) = ty + {1 — t)z, we have for all t £ [0, 1] with the previous estimates 
dPi 



dZjdZq 



{u{t))dt 



< 



i{i-ir\q\'\\Pi\\ C 
Jo 



E 



\ji\>nit)\---\j2nnj,.2it)\dt. 



jeZ^-'^ ,M(j)=~ta~5b 
Multiplying by {zq — yq) and summing in k and j, we obtain 

\\ZpXz) - XpMW,. < Hi - m\P^\\ ( /' \Ht)\g^ dt) \\z - y\\^, . 

J 

Hence we obtain the result after summing in i, using the fact that 
\ty + (1 - 1 < max(||?/|| , ||z|| J. 

Assume now that P and Q are homogeneous polynomials of degrees k and I 
respectively and with coefficients Ofc, h E and 6^, £ E X^. It is clear that 
{P, Q} is a monomial of degree A; + £ — 2 satisfying the zero momentum condition. 
Furthermore writing 

{P,Q}(z)= ''i^i^ 

Cj expresses as a sum of coefficients a^bi for which there exists an a E and 
e E {±1} such that 

(a, e) C fc E Xfc and (a, — e) C I £ Ii, 

and such that if for instance (a,e) = /ci and (a, — e) = £i, we necessarily have 
(A;2, . . . ,kk,i2, ■ ■ ■ ,^1) = J- Hence for a given ji, the zero momentum condition 
on k and on £ determines the value of ea which in turn determines the value of 
(e, a) when = N*^ and determines two possible value of (e, a) when Af = 1,^ . 
This proves (2.11) for monomials. If 

k I 

P = Yp^ and Q = Y Qj 

i=2 j=2 

where Pi and Qj are homogeneous polynomials of degree i and j respectively, 
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then we have 

n=2 i+'j — 2=n 

Hence by definition of ||P|| (see (2.7)) and the fact that all the polynomials 
{Pj, Qj} in the sum are homogeneous or degree i + j — 2, we have by the previous 
calculations 

k+l~2 k+e-2 

\\p\\= E II E {PuQj}\\<'^ E E ^^iii^^ii iiQiii 

n=2 i+j—2=n n=2 i+j—2=n 

k I 

<2fc£(^||p,||)(^||Q,||) =2fc^||p|| iigii 

i=2 j=2 

where we used (2.12) for the last equality. 

The last assertion, as well as the fact that the Poisson bracket of two real Hamil- 
tonian is real, follow immediately from the definition of the Poisson bracket. ■ 



Remark 2.10 The zero momentum condition was crucially used to prove (2.9) 
and (2.11). For instance instead of (2.15), we would have without this condition^ 
in the case where s = 0, 

and this last expression cannot be controlled by the £^ norm of z without an extra 
decreasing property on the \aji\. For instance we can assume that l^'^l 
uniformly bounded with respect to j. See also [15] in the case where the phase 
space is l1 instead of l\ . 

With the previous notations, we consider in the following Hamiltonian func- 
tions of the form 

H{Z) = Hoiz) + P{Z) = ^ LOalaiz) + P(z), (2.16) 

aeAf 

where for all a £ M, 

Ia{z) = CaVa 

are the actions and G R are the associated frequencies. We assume 

VaGAA, \uJa\<C\ar (2.17) 
^Here we assumed that \aj \ does not depend on permutation on the multi index j. 
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for some constants C > 

L 

< 

2.2 Examples 

In this subsection, we present two examples of equations that can be put under 
the previous form. We mention that many other systems can be put under the 
previous form. In particular, we stress out that we do not need any non resonance 
assumption on the frequencies iVa as required in an perturbative approach (see 
for instance [15, 4]). 

2.2.1 Nonlinear Schrodinger equation 

We first consider non linear Schrodinger equations of the form 

idttp = -A^lJ + ^2gi^p,iJ), xGT'^ (2.19) 

where 5 : ^ C is a polynomial of order rg. We assume that g{z,z) G M, and 
that g{z,z) = Odzp). The corresponding Hamiltonian functional is given by 

H{^p,i;) = [ (|VVj|2 + 5(V^,V^)) dx. 

Let (j)a{x) = e*"'^, a G Z'' be the Fourier basis on L^(T'^). With the notation 

/ I \d/2 _ / 1 \ d/2 

^ = [2^) 1^ Ca4>a{x) and V = ) 2^ rjaMx) , 

the Hamiltonian associated with the equation (2.19) can (formally) be written 

ro 

Hit V)=Y1 + 5^ 5Z ^ab Cai • • • ia,m, " " " % ■ (2-20) 

aeZ'* ''=3 a,b 

Here u)a = satisfying (2.17) with m = 2, are the eigenvalues of the Laplace 
operator —A. As previously seen in Examples 2.1, 2.2, 2.7, the nonlinearity P = 
Jfd g{i^{x), 4){x))(lx is real, satisfies the zero momentum condition and belongs to 
Tig (as g is polynomial). 

In this situation, working in the space for ^ corresponds to working in a 
subspace of bounded functions il){x). Similarly the control of the l\ norm of ^ for 
s > leads to a control of HV^Vll^^oo- 



and m > 0. The Hamiltonian system (2.3) then reads 
dP 

Qj^" (2.18) 
= iujaVa + i-^itv), a£j\f. 

eta 
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2.2.2 Nonlinear wave equation 

As a second concrete example we consider a 1-d nonlinear wave equation 

utt - u^x = g{u) , X E (0, vr) , t E R , (2.21) 

with Dirichlet boundary condition: u{0, t) = u{-k, t) = for any t. We assume 
that g : M ^ M is polynomial of order vq — 1 with a zero of order two at u = 0. 
Defining v = ut, (2.21) reads 

\vj \u^x + g{u)J 

Furthermore, let H : H^{0,tt) x L^(0,7r) M be defined by 

H{u, v) = j^^ (^v^ + ]^ul + G(u) j dx (2.22) 

where G such that duG = —g is a polynomial of degree tq, then (2.21) expresses 
as a Hamiltonian system 



d. 



A ^ ( l\ f-u^x + duG 

v) \-i oy V V 

= JVu,vH{u,v) (2.23) 

where ~ (^^^^ qJ represents the symplectic structure and where Vu,v = ^y" 

with V^j and V„ denoting the gradient with respect to u and v respectively. 
Let —Af) be the Laplace operator with Dirichlet boundary conditions. Let A = 
(— A/5)^/^. We introduce the variables {p,q) given by 



q := A^l'^u and p := yl"^/^ 



V. 



Then, on i7''(0,7r) x i?''(0,7r) with s > 1/2, the Hamiltonian (2.22) takes the 
form Hq + P with 

Ho{q,p) = ^{{Ap,p)l2 + {Aq,q)L2) (2.24) 

and 

P(g,p) = [ G{A-^/^q)dx. (2.25) 

In this context M = N \ {0}, uja = a, a £ M are the eigenvalues of A and 
4>a = sin ax, a E AA, the associated eigenfunctions. 
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Plugging the decompositions 



li^) = qa4>a{x) and p{x) = PgCpajx) 
into the Hamiltonian functional, we see that it takes the form 

where P is a function of the variables qa- Using the complex coordinates 

= + ^Pa) and ria = -^{Qa - iPa) 

the Hamiltonian function can be written under the form (2.20) with a nonlin- 
earity depending on G. In this case, the space for z = rj) corresponds to 
functions u{x) such that the Fourier transform u{a) = ^ u{x) sm{ax)dx sat- 
isfies {au)aej\f S i^{M) and {a~^u)aej\f S i^{J\f). This implies in particular 
a control of u{x) and dxu{x) in L°°(0, vr). More generally, with a z in some l]. 
space, s G N, is associated a function u{x) such that d^u{x) S L°°(0,7r) for 
k = 0,...,s + l. 



2.3 Splitting schemes 

In this subsection we describe some Splitting schemes to which we will apply 
our technic in the next sections. The standard Lie- Trotter splitting methods for 
the PDE associated with the Hamiltonian Hq + P consists in replacing the flow 
generated by H during the time h (the small time step) by the composition of 
the flows generated by Hq and P during the same time, namely 

'^Ho°^P and $^o$^^. 

As explained in the introduction, it turns out that it is convenient to consider 
more general splitting methods including in particular a regularization in the 
high modes of the linear part. Thus we replace the operator hHo by a more 
general Hamiltonian ^o- Precisely let ah{x) be a real function, depending on 
the stepsize h, satisfying ah{0) = and ah{x) ~ x for small x. We define the 
diagonal operator Aq by the relation 

Vj = (a,(5) G 2^, AoZj = 6ah{hLJa)zj. (2.26) 

For a G M, we set Xa = ah{huJa)- We consider the splitting methods 

^poi^\ and ^\o^% (2.27) 
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where $p is the exact flow associated with the Hamiltonian P, and where is 
defined by the relation 

Vi = (a, 6) G Z, ($Ao(^))j = exp{iS\a)zj. 
We will mainly consider the cases listed in the table 1 below. 



Method 


ahix) 


Splitting 


ah{x) = X 


Splitting + CFL 


ah{x) = xl;,<c{x) 


Mid-split 


ah{x) = 2arctan(x/2) 


Mid-split + CFL 


ah{x) = 2 arctan(x/2)l^<c(2^) 


New scheme (I) 


= /i^ arctan(/i~'^a;) 


New scheme (11) 


^ ^ x + x'^/h'^ 
"^^^^^ ~ l+x/hP + x'^/h'^P 



Table 1: Splitting schemes 



Let us comment these choices. The "mid-split" cases correspond to the ap- 
proximation of the system 

f]a = iuJaVa a G AT, 

by the midpoint rule (see also [1, 26]). Starting from a given point (Ca^^a) 
have by definition for the first equation 

ea = Cl + fhljl y^ = exp(-2.arctan(/.^j2))e0 
which is the solution at time 1 of the system 

= — i2 arctan(/it(Ja/2)^a a£Af, 
fja = i2 arctan(/iu;a/2)?7a a G AA. 
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Thus in this case the Hamiltonian is given by 

^oiC, ^) = X] ^ arctan(/ia;a/2) (.aVa- 

Note that using the relation 

VyGM, I arctan(y) - y| < (2.28) 

we obtain for all a £ M , 

I exp(— z/iWq) — exp(— 2i arctan(/iLiJa/2))| < Ch^uo^ 
for some constant C independent of a. Using the bound (2.17), we get for all z, 

\\<,{z)-^\^{z)\\^,<Ch'\\z\\^,^. 
More generally, we have the following approximation result: 
Lemma 2.11 Assume that the function ah{x) satisfies: 

Vx > 0, \ah{x) -x\< Ch-^x^ (2.29) 
for some constants C > 0, a > and 7 > 2. Then we have 

II^Ho W - ^Aom^, < Ch-'-'^WzW^^,^ . (2.30) 
Proof. For a given a £ M, we have 

\ah{huj) -huja\< ChT'-'^ujl. 
Hence owing to the fact that |e*^ — e*^| < |x — y\ for real x and y, we obtain 

I exp(-z/ia;a) - eyiY>{-iah{huJa)\ < Ch^'^'io'l, 
and this yields the result. ■ 

Commentary 2.12 Under the assumption that P acts on sufficiently high index 
Soholev spaces H'^ , the previous result can he combined with standard convergence 
analysis to show that the splitting methods (2.27) yield consistent approximation 
of the exact solution <I>^ provided the initial solution is sm.ooth enough (depending 
on m). The condition 7 — cr > 2 guarantees a local order 2 in (2.30) that will 
he of the same order of the error made hy the splitting decomposition after one 
step. Such a local error propagates to a glohal error of order 1, which means that 
for a give finite time T , the error with the exact solution after n iterations with 
nh = T will he of order n x c^^ h up to constants depending on T, and under 
the assumption that the numerical solution remains smooth. Give more precise 
results would he out of the scope of this paper, and we refer to [22] for the case of 
NLS. 
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Let us consider the function 



ah(x) = /i^arctan(/i-^x) (2.31) 

for 1 > /3 > 0. It satisfies (2.29) with C = 1/3, a = 2/3 and 7 = 3 (see (2.28)). 
Hence for /? = 1/2, the estimate (2.30) shows a local error of order 7 — cr = 2, 
and hence the splitting schemes (2.27) remains of local order 2 (though with 
more smoothness required than with the midpoint approximation) which means 
that the error made after one step is of order h? . Of course, when (3 = 1, the 
approximation is not consistent (local error of order 1). 
The second example 

cthix) = ^;tr ^Kno» (2-32) 

exhibits similar properties. Note that the simple choice 

a,{x) = ^-p^ (2.33) 

ensures only a local error of order h^^f^ (7 = 2 and a = ;9 in (2.30)). Hence the 
corresponding splitting schemes (2.27) are of global order 1 — (3 (hence 1/2 in the 
case where /? = 1/2). 

All these "new" schemes have the particularity that ah{x) ~ x when x is 
small, but when x ^ 00, we have ah{x) ~ . Their use thus leads to a stronger 
regularization effect in the high modes than the midpoint approximation, with- 
out breaking the order of approximation for smooth functions. We will see in 
section 3 and 4 that this property allows to construct a modified equation over 
exponentially long time for all these schemes. 

Note that in practice, the implementation of the schemes associated with 
the filter functions (2.31) or (2.32) a priori requires the knowledge of a spectral 
decomposition of Hq. This will be the case for NLS on the torus, or NLW with 
Dirichlet boundary conditions, where the switch from the x-space (to calculate 
$p) to the Fourier space (to calculate ^\^) can be easily implemented using fast 
Fourier transformations. 



3 Recursive equations 

In this section we explain the strategy in order to prove the existence of a modified 
energy. We will see that it leads to solve by induction a sort of homological 
equation in the spirit of normal form theory (see for instance [15]). For simplicity, 
we consider only the splitting method <I>pO$^^^ . The second Lie splitting <1>^^ 0$^ 
can be treated similarly. 
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We look for a real Hamiltonian function Z{t, ^, rj) such that for all t < h we 
have 

and such that Z{0) = Aq. 

For a given Hamiltonian K £ TCs, we denote by Ck is the Lie differential 
operator associated with the Hamiltonian vector field X^- for a given function g 
acting on £l, s > 0, and taking values on C or £], we have 



dz 



■J 



Denoting by z{t) the flow generated by Xk starting from z £ l\, i.e. z{t) = $^(z), 
we have (if K is in TLs) 

z^^\t) = C'k[I]{z{t)) for all A; G N, 

where / define the identity vector field: I{z)j = Zj. Thus we can write at least 
formally 

$), = exp(£K)[/]. (3.2) 
Differentiating the exponential map we calculate as in [18, Section HI. 4.1] 

where the differential operator associated with Q{t) is given by 

fc>0 ^ ^' 

with 

AdcA'^H) = [Ca,Ch] 

the commutator of two vector fields. 

As the vector fields are Hamiltonian, we have 

i^A, Ch] = ^{A,H}, 

where 

adx(G) = {K,G}. 
Hence we obtain the formal series equation for Q: 

^(^) = E7r^a4w^'W (3-3) 

fe>0 ^ ^' 
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where Z'{t) denotes the derivative with respect to t of the Hamiltonian function 
Z{t). 

Therefore taking the derivative of (3.1), we obtain 

Xp O $*p O = O 

and hence the equation to be satisfied by Z{t) reads: 

T.jkTTy,^^m^'('^ = P- (3-4) 

Notice that the series J2k>o {k+iy. '^^ ~ ^"T^ invertible in the open disc \z\ < 
2tt with inverse given by X^fc>Q where Bf^ are the Bernoulh numbers. So 
formally, Equation (3.4) is equivalent to the formal series equation (see also [8], 
Eqn. (3.1)) 

^'W=E|f^d|(,)P. (3.5) 

Plugging an Ansatz expansion Z{t) = J2e>o ^^Zi into this equation, we get 
Zq = Aq and for n > 

{n + l)Zn+i = ^^ adzt^ ■ ■ ■ adz.^P. (3.6) 

fc>0 ■ ti-i |-4=n 

Commentary 3.1 The analysis made to obtain this recursive equation is formal. 
To obtain our main result, we will verify that the series we manipulate are in fact 
convergent series in Hs uniformly on balls of £l that contains the different flows 
involved in the formulas (see in particular Lemma 4-3 below). 
For instance (3.2) holds true as soon as z{t) remains in a ball Bfj- := {z S 

■^s I IklLi ^ M} for < t < 1, and that the series ^^L?^^"* ^•^ uniformly 
convergent on B'^j. Notice that this in turn implies that t z{t) is analytic on 
the complex disc of radius 1. 

Commentary 3.2 In the case of Strang splitting methods of the form 

we can apply the same strategy and look for a Hamiltonian Z{t) satisfying for all 
t < h, 

and this yields to equations similar to (3.6). We do not give the details here, as 
the analysis will be very similar. 
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The key lemma in order to prove that the previous series converge (and thus 
to justify the previous formal analysis) is the following one (and whose proof is 
straightforward calculus) : 

Lemma 3.3 Assume that 

Q{z) = ^ ajZj 

is a polynomial, then 

adzo(Q) = ^ik{3)ajZj 

where for a multi-index j = {ji, ■ ■ ■ , jr) with for i = l,...,r, ji = {ai,6i) G 
J\f X {±1}, we set 

Hence we see that if A(j) < 27r we will be able to define at least the first term 
Zi by summing the series in k in the formula (3.6). 



4 Analytic estimates 

We assume in this section that Oh and h satisfy the following condition: there 
exist r and a constant 6 < 2tt such that 

Vn<r yjein, |A(j)| < 27r-5. (4.1) 

Using Lemma 3.3, this condition implies that for any polynomial Q G Vr, we 
have the estimate 

lladzoQII < (27r-5)||Q|| (4.2) 

as for homogeneous polynomials, the degree of adz^Q is the same as the degree 
of Q. 

Theorem 4.1 Let ro > 3. Assume that P £ Pro and that the condition (4.1) 
is fulfilled for some constants S and r. Then for n < N := we can define 

polynomials Z„ G Vn(ro-2)+2 satisfying the equations (3.6) up to the order n, and 
satisfying the estimates \\Zi\\ < c and for 2 < n < N , 

\\Zn\\ < c{CnY"'^ (4.3) 
for some constants c and C depending only on \\P\\ , ro and 5. If moreover 

P G SPro then Zn G SPn{r^_2)+2- 
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Proof. Let 

1=1 jail. 

First we prove the existence of the for k < N. The equation (3.6) for n = 
reads 

k>0 

The previous Lemma and the condition (4.1) show that Zi exits and is given by 

Further we deduce immediately that Zi is real, and satisfies ||.Z'i|| < qHPH for 
some constant q. 

Assume now that the are constructed forO<A;<n,n>l and are such that 
Zk is a polynomial of degree k{rQ — 2) + 2. Formally Z^+i is defined as a series 



n + - 

A;>0 

where 

Ak= ^ ad^^^ • • • adzf.^ P. 

ei+-+ik=n 

Let us prove that this series converges absolutely. In the previous sum, we sepa- 
rate the number of indices j for which ij = 0. For them, we can use (4.2). Only 
for the other indices, we will use the estimates of Proposition 2.8 by taking into 
account that the right-hand side is a sum of terms that are all real polynomials 

of degree (4 H h 4)(ro - 2) + ro = {n + l){ro - 2) + 2 < {n + l)ro and hence 

the inequality of Proposition 2.8 is only used with polynomials of order less than 
(n -|- l)rQ. Thus we write for A; > n 

Pfcll := II Yl ^^^^1 •••adz.^^ll < 



i=i ' ' ei+-+ei=n\ej>o 

n 

< {27T- 6)'-^ k^Yl E (n + l)^-i2V2^4||%|| •••^ill^^ll ||i= 
i=i e^j^...j^e.=n\£j>o 

and thus X]fc>o ^^^^ converges and is well defined up to n -|- 1 < A^. 
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Further, as the right-hand side of (3.6) is a sum of term that are all real polynomi- 
als of degree (n + l)(ro — 2) + 2, Zn+i is a polynomial of degree (n + l)(ro — 2) -f2. 
Now we have to prove the estimate (4.2). Following the previous calculation we 
get 

(n+l)||Z„+i|| < 

n „ 

1=1 £iH \-ei=n\ej>0 k>i ^ ' 

h\\Zi,\\ ■ ■ ■ iiWZi^W \\p\\ . 

On the other hand, the entire series f{z) := J2k>i 1^^^ defines an analytic 
function on the disc \z\ < 27r. Thus its i — th derivative X]A:>i (^k-i)\ ^^~^ 
defines an analytic function on the same disc and, by Cauchy estimates, there 
exist a constant Cs = sup|2|<27r_5/2 1/(2^)1 such that 

We then define for n > 1 

Cn = nro||Z„|| . 

These numbers satisfy the estimates, for n > 0, 

n 

Cn+1 < + C^rollPlI E + l)'"'(4rVo)X,, • • • C^. 

i=l ^j+...+4=n|£j>0 

Let us fix > 1. We have for n = 0, . . . , 

4rVo(iv + i)c„+i < 4rVo(iv + i)6%cs\\p\\ 

n 

+ csA5-\i\\p\\ J2 E + ir(4rVo)x^i ■■■Ci.. 

i=l £j+...+£.=„|£^.>o 

Let Pj, j = 0, . . . ,N the sequence satisfying 

n 

1 = 1 £j + ...+£.=„|^^.>0 

where 

Ci = 4(iV + l)r2c5 r^PlI and C2 = ^r2||P|| . 
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By induction, we have that for all n = 0, . . . , A^, 



(iV + l)^Cn,</3n. 

Multiplying (4.4) by t"-+^ and summing in n > we see that the formal series 
f3{t) = Ylj>i ^"'/^i satisfies the relation 

This yields 

{i-p{t)){m-tci) = tc2P{t) 

or equivalently 

- p{t){l + t{Ci - C2)) + tCi = 0. 
The discriminant of this equation is 

(1 + t{Ci - C2)f - 4tC7i = 1 - 2t(Ci + C2) + t\Ci - C2f 
and hence, for t < l/2(Ci + C2), we find using /3(0) = 0, 

2/3 = 1 + t{Ci - C2) - (1 - 2t(Ci + C2) + t\Ci - 02)"^ ^''^ 



We verify that for t < l/2(Ci - C2) we have 

21} < I. 

By analytic estimate, we obtain that for all n > we have 

n! 2 V 
For n = N , this yields 

l5N<{CN\\P\\f 

for some constant C depending on ro and 5. We deduce the claimed result from 
the expression of Cat. ■ 

For s > 0, we define 

Blt = {z^eMz\\^,<M]. 
and we will use the notation Bm = B^. 

Theorem 4.2 Let tq > 3, s > and M > I be fixed. We assume that P £ Pro 
and that the condition (4.1) is fidfiUed for some constants 6 and r > and we 
denote by N the largest integer smaller than ^~_^2 ' ^^^n there exist constants cq 
and C depending on tq, s, 6, \\P\\ and M such that for all hN < cq, there exists 
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a real Hamiltonian polynomial G VN{rQ-2)+2 such that for all z G we 
have 

\\^% o ^\^{z) - < h'^+^CNf. (4.5) 

Moreover, assuming that 

ro 

£=1 jeli 

then for z G B^^ we have 

\Hhiz) - Hj^\z)\ < Ch (4.6) 

where 

aGAT e=i jell " 

If finally, P G SVro, ^^^^ G '5Pyv(ro-2)+2- 
Proof. We define the real Hamiltonian Hh = — 7^-^, where 

AT 

j=0 

and where, for j = 0, ■ ■ ■ , N, the polynomials Zj are defined in Theorem 4.1. 
Notice that A'' depends on r and thus on h via the condition (4.1). 
By definition, Zh{t){z) is a polynomial of order N{rQ — 2) + 2 < NrQ and using 
Theorem 4.1 we get 

N 

\\Zhit)\\ <ci(l + ^(Ctj7-2)<oo. 
J =2 

for some constant ci depending on 6 and Hi-"!!. Thus Z^ G 7^jv(rQ_2)+2 (and in 
5P^(,0_2)+2 if PGcsnj. 

We will use the fact that for all s, there exists a constant Cg such that for all 
J >1, 

f < {csy. (4.8) 

Now, as for all j, Zj is a polynomial of order j{ro — 2) + 2 < jVo with a zero 
of order at least 2 in the origin, we have using Proposition 2.8 and Theorem 4.1 
that for z G Bgj^j and j > 1. 

||Xz,(z)||,, < 2c{jror^\Cjy-'{ sup ||z||' ) < M{2Ccy,-^'cj{9MYy 

fc=l,...,jro-l 
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where the constants c and C are given by estimate (4.3). On the other hand we 
have using Lemma 3.3 and (4.1) 

||Xzo(^)ll,i <27r||z||^, <27rM. 
Hence, for t < (4iVCc,r^+^c)^i(9Af)-'^o we have 

N 

\\Xz^^t)iz)\\^, < 27rM + M^{t2Ccsr'Q-^^cN{9MY'>y < (27r + l)M < 8M. (4.9) 

i=i 

Therefore by a classical bootstrap argument, the time 1 flow ^^^(t) ^ti 

provided that t < (4ArCc,r^+^c)-H9M)-^o. 
On the other hand, is an isometry of and hence maps B^^ into itself, 
while using again Proposition 2.8, we see that $p maps B'^j into Bgj^ as long as 
t < (4ro+^||P|| M(^o^i))-i. We then define 

T^T{N,M,ro,s,6,\\P\\) 

:= min{(4r^+^||P|| M^''°-^Y\ (4iVC7c,r^+^c)-^(9M)-'"°} (4.10) 

and we assume in the sequel that < t < T in such a way that all the flows 
remain in the ball Bqm- 

Let u{t) = o ^^^(z) — <I>^ denote by Qhit) the Hamiltonian defined 

by 

fc>0 ^ ^' 



Lemma 4.3 Fort < T given in (4.10), the Hamiltonian Qh{t) G 'Hs and satisfies 
for zeBlj 

^^kw(-)=^Q.W°^kw(^)- (4-11) 

We postpone the proof of this Lemma to the end of this section. 
Using this result, as u(0) = 0, we get for t < T given in (4.10) 

11^(011,, < /"Vp('J>po<^Ao(^))-^Q.{s)(^k{s)(^))l^ids 


< f ||Xp($^^(,)(z)) - XQ,^(.)(a>^^(,)(z))||,, d. 
+ f \\Xp{^% o ^\^{z)) - Xp(^.^^(,)(z))||^, ds. 
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Therefore for t < T 



\u 



„< f sup \\Xp{z)-XQ^i^,-^{z)\\^,ds + Lp f \\u{s)\\ ds (4.12) 

Jo ZdBgM " Jo 

where using Equation (2.10) in Proposition 2.8, we can take 

Lp = 4rg+i||P|| (9M)'"o-2. 

So it remains to estimate sup^go ||Xp(z) — Xq^(j)(^;)|L^ for z E Bqm and t <T. 

3 

Now by definition of Qh{t) and using Lemma 4.3 we have 

oo „ 

fc=0 

where the right hand side actually defines a convergent series by the argument 
used in the proof of Theorem 4.1. By construction (cf. section 3), we have 

oo „ 

k=0 

in the sense of Hamiltonian in the space 7is- Taking the inverse of the series, we 
see 

Qhit) -P=Y,Kn (4.13) 

n>N 

where we have the explicit expressions 

K. = (m+1)^ Yl adz,^ •••adz,^Z^+i. 



(A; + _ 

l+m=n\ m<N k>0 ^ '' ^iH h4=^Kj<JV 

Estimates similar to the one in the proof of Theorem 4.1 lead to 

\\^n\\< E (^+i)E2vg>+i)-E ^';r/_jy >< 

l+m=n\m<N i=0 k>i 

h\\Zi-i\\ • • • ii-i\\Zi._-^\\ £i\\Z£.\\ \\Zm+i\ 

£i+...+£i=£|0<£j<Af 

and hence after summing in fc, 



(4.14) 



\Kn\\<Ci E ("^ + 1)E 

i+m=n\ m<N i=0 



2Vf(n + ir 

: 



E h\\Zi-^\\ ■ ■ ■ £i-i\\Zi._-^\\ ii\\Ze-\\ \\Z., 

£l+...+£i=£|0<£j<Af 



m+l I 
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for some constant Ci depending on 6. Using the estimates in Theorem 4.1, we 



have for £j > and ij < N, 



Using moreover \\Zm+i \\ < c^CN)"^^ and as the number of integer £i, . . . , £i stricly 
positive such that + = £ is bounded by 2^^, we obtain 



\Kn\\<cCi 



(2c)V^*(?i + 1)* 



l+m=n m<N i=0 



ei+-+ii=e\o<e^<N 



£+m=n i=0 

Therefore, there exist a constant D depending on ro and ||P|| such that 

Vn > iV, ||K„|| < (DNY. 

As ii'n is a polynomial of order at most ron, we deduce from the previous estimate 
and Proposition 2.8 that, for z G B^j^^ 

\\XkA^)\\^, < linroY^HDNrmr^ < {2csr'+^DN{9MY°r, 

where the constant Cg is defined in (4.8). Hence the series X^.„>o ^^-^^XnC-z) con- 
verges for t < {AcsrQ-^^DNidMY'o)'^. Furthermore, again for z E and 
t < {4.Csr'Q+^DN{9MY")~^, we get using (4.13) and the previous bound 

||Xq,(,)(z)-Xp(z)||^, <Y,nXKAm,^<{N + l)t''{BNf 

n>N 

for some constant B depending on ||P|| , s, 6 and tq. 
Let us set 

Co{M,ro,6,s,\\P\\) : = 

min{(4r^+i|||P|| M'-''°^'^^)~\ (4Cc,r^+^c)-i(9M)"^o , (4c,r^+^Z)(9M)''«)-i}. 
For t < co{M,rQ, 6, s, ||P|| )A^~-'^, inserting the last estimate in (4.12) we get 



and this leads to 



<t^+\BNf + Lp ['\\n{s)l,ds. 
Jo 
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for some constant B depending on ro, (5, s \\P\\ and M. This implies (4.5) defining 

Hh = Zh{h)/h for h < co{M,ro,S,s, ||P||)A^-^ 

The second assertion of the theorem is just a calculus defining 

Using the previous bounds and the first inequality in Proposition 2.8, we then 
calculate that for z G B^^j 

N 

N 

i=2 
i=2 

by definition of cq{M, vq, 6, s, \\P\\ ). ■ 
Proof of Lemma 4.3. With the previous notations, we have 

n>0 

where Kn is given by (4.14) and the bounds obtained show that Qh(t){z) and 
-^Qh{t)(^) ^^11 defined on Bqm- Now let us consider the flow '^'^/(t)- 
previously mentioned, it acts from Bf^ to Bg^^. Now we can write formally 

k>0 



n>0 k>0 ei+-+ek=n\ei<N 



(4.15) 



n>0 



We are going to show that this series converges uniformly for z £ B\j and t <T . 
Let K be fixed polynomial of degree /c, and G{z) = {Gj{z))j^z a function acting 
on £l, and taking value in il, and such that the entries Gj{z) are all polynomials 
of degree i. By definition, we have 

8C 

{Ck o G), = Y,{Xk\-q^ = {K, G,}. 
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Now using the relation (2.11) of Proposition 2.8, we have that {Ck ° G)j is a 
polynomial of degree k + i — 2 and of norm smaller that 2A;£||i('|| . Now if 
K = Zq, this bounds can be refined in \\{Czo o G)j\\ < (2'/r — || using (4.2). 

For a given z £ the j-th component of jCzi_^ o • • • o C^e^ is a polynomial 

of order n(ro — 2) + 2 and involve terms of momentum M.{j) = —ea if j = (a, e) 
(see (2.13)). Hence summing in j, and separating as before the indices m for 
which £m = in the sum, we obtain for a given n and z G 

i=l k>iii+---+ii=n\en<N ^ ' 

and we conclude as before that this series is convergent for t < T given in (4.10) 
with K depending on tq, M, ||P|| , s and 5. 

Now writing down the same argument for the series (in A; > and t", n > 0) 
defining ^^^^(4) and Xq^(i-^ thd^t this series are again convergent, 

which justify the relation (4.11). ■ 



5 Applications 

In this section we analyze the consequences of the analytic estimates obtained 
in the previous section. We first show that for the "new schemes" in table 1, we 
obtain exponential estimates. We then show that for general splitting schemes, 
we obtain results under an additional CFL condition depending on r. 

5.1 Exponential estimates 

We consider the following splitting scheme 

^pO^Ao (5.1) 
where the operator Aq is associated with a function ah{x) (see (2.26)) satisfying 

VxGR |ah(a;)| < 7/1^ (5.2) 

for some /? G (0, 1) and some 7 > 0. Examples of such methods preserving the 
global order 1 approximation property of the scheme (5.1) for smooth functions 
are given in table 1. For such scheme, we obtain an exponentially closed modified 
energy: 

Theorem 5.1 Let ro > 3 s > and M > 1 be fixed. Assume that P £ Vtq 
and that ah satisfies the condition (5.2) for some constants 7 and j3. Then there 
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exists a constant ho depending on ro, ||P|| , s, M and 7 such that for all h < ho, 
there exists a real polynomial Hamiltonian Hh such that for all z S we have 

\\^% o ^\^{z) - ^^uSm,. < hexp{-{ho/hf). (5.3) 

Proof. The hypothesis (5.2) imphes that the eigenvalues Xa of the operator Ao 
are bounded by 7/1^. Hence for a multi-index j = (ji, . . . ,jr) we have 

|A(j)| < r7/i^ 

and the condition |A(jf)| < vr will be satisfied as long as r < (7r/7)/i-^. Taking 
ri such that ri < {7r/^)h~^ < ri + 1 and defining N = {ri — l)/(ro — 1), we get 
bih~^ < N < b2h~^ for some positive constants 61 and 62 depending on 7 and 
ro- Now the estimate (4.5) in Theorem 4.2 for this yields 

11$^ o ^\^^{z) - ^'H^m,^ < h'^+HCNf < h{Cb2h'-T, 

as long as hN < cq. Thus defining 

ho = min{{cob^y/('-^\{eCb,)-'/^'-^\bl/^}, 
we have for < /i < /iq 

||<i>^ o <^>X^{z) - ^'h.mei < he-"" < he-"^^-' < heM-iho/hf). 

■ 

The dynamical consequences for the associated numerical scheme are given in 
the following corollaries. We first assume that the numerical solution remains a 
priori in only over arbitrary long time. 

Corollary 5.2 Under the hypothesis of the previous Theorem, let z^ = (C°,C°) G 
and the sequence z^ defined by 

= $^o$\^(z"), n>0. (5.4) 

Assume that for all n, the numerical solution z'^ remains in a hall Bm of for a 
given M > 0. Then there exist constants ho and c such that for all h < ho, there 
exists a polynomial Hamiltonian such that 

Hh{z^) = Hh{z'^) + 0{eM-ch-^)) 

for nh < exp{ch~^). Moreover, with the Hamiltonian H^^ defined in (4.7) we 
have 

H^\z-) = H\l\z^) + 0{h) (5.5) 
over exponentially long time nh < exp(c/i^^). 
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This means that the modified energy remains exponentially closed from its initial 
value during exponential times. Or more practically (since i?^^^ is explicit) the 
first modified energy is almost conserved during exponential times. 
Proof. As all the Hamiltonian function considered are real, we have for all n, 
_ (^"j^")^ i.e. z"' is real. Hence for all n, Hh{z"-) G M. Note moreover that 
we can always assume that M > 1. 

We use the notations of the previous theorems and we notice that Hh{z) is a 
conserved quantity by the flow generated by Hh. Therefore we have 

Hh{z-+^) - H,{z^^) = H,,{^% o ^\^{z-)) - H,i<P%^izn) 

and hence 

mz-+') - H,{zn\ < ( sup ||Vi7,(z)||^^)||cI>^oci>3^^(z-)-$^J.-)||^,. 

Now using (4.9) and (5.3) we obtain for all n 

\Hh{.z^^^) - Hhiz^l < A7rMhexp{-{ho/hf) 

and hence 

\Hh{z^) - Hhiz'')] < (n/i) exp(-2c/i-^) 

for some constant c, provided /iq is small enough. This implies the result. The 
second estimate is then a clear consequence of (4.6). ■ 

The preservation of the Hamiltonian H^'^ over long time induces that for 

z" = (r,<), 

aeAf 

is bounded over long time, provided z^ is smooth. For the functions ah given 
in Table 1, this yields i/™/^ bounds for low modes, while the norm of high 
modes remain small (see [8, Corollary 2.4] for similar results in the linear case). 
We detail here the result in the specific situation where ah{x) is given by (1.10). 



Corollary 5.3 Let tq > 3, P G Vtq and ah{x) = arctan(x/\//i). We assume 
that there exists a constant 6 > 1 such that 

ya£M, ^\ar < iOa < b\ar . (5.6) 

Let the sequence defined by (5.4). We assume that z^ = {(,^,^^) G H"^/"^ and 
that the sequence z^ remains hounded in i^. Then there exists constants C , a and 
(3 such that 

over exponentially long time nh < exp(c/i~^). 
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Proof. The almost conservation of the Hamiltonian Hj^^ (cf. (5.5) ) shows that 
for all n such that nh < exp(c/i~^), we have 

V ^arctan(^/)Iwa)|eaP 

< \Pt'\z^) - + E ^ arctan(^..,)|e,T + Ch 

for some constant C, where P^\z) is the non linear Hamiltonian of equation 

(4.7). As remains bounded in , we have that \P^\z^) — P^^\z'^)\ is uniformly 
bounded. Now we have for all a 

— p arctan(v7iWa) < ^a- 

and we deduce using (5.6) 

V ^arctan(\//^!wa)|^°|2 < C\\z\\ . 

Hence we have 

V 4farctan(^/^cuj|ei' < C 

for some constant C depending on M and ||2||^„j/2- We conclude by using 
X > 1 =^ arctanx > arctan(l) and x <1 =^ arctanx > — . 

and the bounds (5.6) on Wq. ■ 

Now we conclude this section by considering the case where the numerical 
solution z^ remains a priori bounded in some s > 0. In this case, we can show 
that the initial energy H{z^) remains bounded over exponentially long time, 
provided that s is sufficiently large (typically of order 3 in concrete applications). 
We first begin with the following Lemma: 

Lemma 5.4 Assume that h < 1, and that ah{x) satisfies: 

Vx>0, \ah{x) -x\ <Ch^''x'^, (5.7) 

for some constants a>0,^>2 + a and C > 0. Assume moreover the non reso- 
nance condition (4.1) is satisfied for some 6 > 0, and let H^^^ be the Hamiltonian 
defined in (4.7). Then for all M > and s > (where m is given in (2.17)J, 

and for all z £ Bf^, 

\H;^\z)-H{z)\ <Cih 
where Ci depends on M and 5. 
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Proof. For all a G M, we have using (5.7) and 7 — o" — 1 > 1 



-ah{huJa) - 
h 



Hence as 2; = rj) is real, 



— Hq{z) 



<ChY, ^liaria < Ch\\z\\], < CM^h. 



Now, for all |x| <2'k — 5, we have 



1 



< Qle'^ -l-ix\< Cs\xf, 



e"' - 1 

for some constant Cs depending on 6. Hence we have for all z G il, 



EE 



\Zn\. 



1=1 jele 

Now using the fact that for a given multiindex j = {ji, . . . , jr^ ) with for all 
i = 1, . . . , ro, ji = (oj, 6i) , we have 

ro 

|A(j)| < roh max Iw^J < ro/iTT|iir, 

i=l,...,ro 



1=1 



we obtain 



ro 



EE 



exp(iA(j)) -1 ' ' ^ ' 



< Csi27r - 6)roh\\zC 



As 7717/2 > 77T,, we have \\z\\^° < IvF° , and this yields the result. ■ 

As an example of application, we give the following result: 

Corollary 5.5 Under the hypothesis of Corollary 5.2, assume that ah{x) = 
\//iarctan(x/\//i), and that the sequence defined by (5.4) remains in a ball 
of il for a given M > and s > 3m/2 where m is given in (2.17). Then 
there exist constants Iiq and c such that for all h < ho, and all n such that 
nh < exp(c/\//i), we have 



H{z'^) = H{z°) + 0{h) 
where H = Hq + P is the original Hamiltonian. 



(5.8) 
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Proof. In view of Corollary 5.2, we only have to prove that for real z G 
we have 

\H^^^ {£) - H{z)\ <Ch 
which is given by the previous Lemma with 7 = 3, cr = 1 and 5 = vr. ■ 



5.2 Results under CFL condition 

We now consider cases where an does not depend on h and thus does not satisfy 
(5.2). We focus on schemes such that 

ah{x) = ^(x)l^.<c(x) (5.9) 

i.e. schemes associated with a filter function /3 and a CFL condition with CFL 
number c. We mainly have in mind the two applications (see Table 1) 

/3(x) = X and /3(x) = 2arctan(x/2) 

corresponding to standard and implicit-explicit splitting schemes. Now for such 
scheme and for all multi-index j = {ji, ■ ■ ■ ,jr) we have at most 

\A{j)\<r(3{c). 

Hence if we define 

Cr = r\^), (5.10) 

the condition (4.1) will be satisfied for some 6 if c < Cr. The following result is 
then a easy consequence of Theorem 4.2: 

Theorem 5.6 Let vq > 3, r > ro, s > and M > be fixed. Assume that 
P G "Pro c-i^d that ah is of the form, (5.9) for som,e constants c < Cr where Cr 
is given by the equation (5.10). Then there exist constants ho and Cr such that 
for all h < ho, there exists a real polynomial Hamiltonian such that for all 
z G Bf^, we have 

ii$^o$i,^(z)-ci>^j^)ii^, <a/^^(^) (5.11) 

where N{r) = (r - 2)/(ro - 2) + 1. 

// moreover P G SVr^ and for all a G Af, oJa > 0, then the same results holds 
true under the weaker condition 

C<Cr/2=r\ — )- (5.12) 

r 
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Proof. The first part of the Theorem is a consequence of the previous estimates. 
The last assertion comes from the fact that all the Xa are positive. Hence for a 
given r and a given monomial of SVr associated with a symmetric multi-index 
j = (ji, • • ■,jr/2,h, K/2) with ji = (ui, +1) and h = {bi, -1), and bi £ M, 
we have 

r/3(c) 

r/2 ''"1 • • • '"Or/2 ~" 



= Aai + . . . + Xa,.,2 - Xbi- Xb < 



and this yields the bound (5.12). ■ 

This result implies the preservation of the Hamiltonian over long times of 
order h~^^'^^^^, as in Corollary 5.2. Similarly, Corollary 5.3 extends to the case 
of the mid-split scheme (i.e. the case (3{x) = 2 arctan(a;/2)) leading to a i/"^/^ 
control of the frequencies smaller than Crh~^^"^, over long times of order , 
namely 

\jr\z^\'<c 

for n < h~^^'^\ As the proof is completely similar, we do not give the details 
here. Note however that the high frequency cut-off does not give a control of the 
norm for high modes. 

In Table 2, we give the expression of the CFL constant Cr in (5.10) required in 
order to obtain a given precision of order h^^^^ in the estimate (5.11) in the two 
cases where (3{x) = x and f3{x) = 2 arctan(x/2). The third column is concerned 
with the nonlinear cubic Schrodinger equation with a general cubic nonlinearity 
(i.e. a quartic Hamiltonian, ro = 4) while, in the fourth one, we consider the 
case of the nonlinear cubic Schrodinger equation with a nonlinearity belonging to 
SV4. For instance, we see that in this latter case, the preservation of the energy 
is ensured over long times of order with a CFL or oder 1 for the mid-split 
integrator. 

Commentary 5.7 Note that Corollary 5.5 can be easily extended to the situa- 
tion where (3{x) = 2 arctan(2;/2) with a CFL condition Cr- Indeed, in this case, 
for a real z = (i^,?]) the difference between the energy Hj^\z) with ah{x) = 
f3{x)lx<:cj.{x) and the same energy but associated with ah{x) = /3(x) is 

2arctan(/ia;a/2)^a??a < vr ^ Ca??a 

a&Af\ hua>Cr aSA^I h|a|'">Cr 

2 

and this last expression is bounded by Ch\\z\\ for s = m/2. We then conclude 

using Lemma 5.4 with ah{x) = /3(x), a = and^ = 3 to prove statements similar 
to Lemma 5.5 in this situation. 
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B(t] = T 


B(t] = 2 arrtan I T /2 1 


(3{x) = 2 arctan(a;/2) 
for cubic NLS 




1.57 


2.00 


oo 




1.05 


1.15 


3.46 




0.79 


0.83 


2.00 




0.63 


0.65 


1.45 




0.52 


0.54 


1.15 




0.45 


0.46 


0.96 




0.39 


0.40 


0.83 




0.35 


0.35 


0.73 




0.31 


0.32 


0.65 



Table 2: CFL conditions for quartic nonlinearities 



6 Equations on R 



d 



In this section we generalize our method to the "continuous" case, i.e. the case 
where the hnear part of the Hamiltonian has a continuous spectrum. As an 
example we consider the nonlinear Schrodinger equation (NLS) on R*^: 

iut = —A.U + d2g{u,u) (6.1) 

where we assume that 5 : — > C is smooth and has a zero of order at least 3 at 
the origin and satisfies g{z, z) S M. 

The natural phase spaces are the Sobolev spaces, {u,u) G H'^{W^) x H'^iW^) 
which are endow with the symplectic structure associated with the following 
Poisson bracket 

{FiG}{u,u) i^du du du du) 

In this Hamiltonian setting the NLS equation reads 



dH _ dH 

lUt = 

with 



ou ou 



H{u,u) = / [\'Vu{x)\'^ + g{u{x),u{x))) dx. 

The Fourier transform maps the symplectic space H'{R'^) X H'{R'^) into the 
symplectic space := L^(M'^) x LKR"^) where 

L'M') ■■= W e / (1 + \C\'im)\'dC < +00} 
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and, in the Fourier variables, the Poisson bracket reads 

In this setting the NLS reads 

i4>t = 



dip dip 



dP 



OH 

dip 

dH _ ^2 dP 



dip' 



with 



and 



H{<p,ip)= f emmdx+p{<p,^) 



mV')=/ gi 0(e)e""-«de, / V(??)e^"-''dr? ) dx. 
In particular for the cubic NLS, g{u,u) = and 



P{(P,^)= / '/>(6)'/'fe)V'(^i)^(%)d6de2dr/id772. 

>^?i+?2-r?l-»?2=0 

More generally, i.e. for general (analytic) nonlinearity 5, 

P{<P,iP) = 



where 
and 



yi "i/ / ■ ■ ■ 4>{ij)'^{'>n) ■ ■ ■ ^(%)d6 • • • d^jdr/i • • • di]i 

??) = 6 + • • • + - m m 



To give a general framework similar to the discrete case we introduce the following 
notations: for a given d, we set 



x{±l} 



endowed with the canonical measure d^ x For r = {£,-,5) G Z'^^ we set 

|r| = |.^|, and we define the variable z{t) = ((/>(^), V'(C)) by the formula 

z{i,5) = m if 5 = 1, 

z{i,6) = iPiO if <5 = -l. 
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For s > 0, we say that z belongs to L\ if the norm 

\\A\^. = \ [ a + mi\m\ + \m\)dc 

is finite. 

Following Definition 2.3, we define 

Definition 6.1 We denote by TC^ the space of Hamiltonians P satisfying 
PeC^{L\X), and Xp GC^iLlLl). 

For instance 

Piz) = [ </'(ei)</'(e2)V'(r/i)V'(%)deide2dr?idr?2 
is in Tig with 

Xp{z)it 1) = 2i / m)mMv)d^id^2dv 

and 

Xp{z){C, -1) = -2i I </'(ei)V'(??i)V'(%)dadr?id772. 

As in the discrete case, we say that z is real if z = (</), and that a Hamiltonian 
P{z) is real if it is real for real z (like the example just above). 

For the multi- variable r = (ti, . . . , r^) G {Z'^Y we use the notation 

z{t) = z{ti) ■ ■ ■ z{Ti). 

We define the momentum M{t) of the multi- variable r by 

M{t) = ii5i + ■ ■ ■ + itS^. 

We then define the subset of {Z'^Y 

lf = {T = {Ti,...,n)(^Z\ with M{t) = {)}. 

We also use the notation r = {$,,6) with ^ = (^i,...,^^) and 5 = [5i,...,5ii) 
which means that for all i = 1, . . . , ^ we have ti = (^j, 5i). We denote by dlf{r) 
the measure defined by 

where d^ = d^i • • • d^^. 
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Definition 6.2 We say that a polynomial Hamiltonian P € if it is real and 
if we can write 



with 



P{z) = j2[ a{T)z{T)dJt{T) 
p='). •''^i 



El|o|Uoo(x^=: ll^ll <oc. 



Similarly to the discrete case, we define the set of symmetric polynomials SVk- 
The following result is the analog of Proposition 2.8: 

Proposition 6.3 Let k >2 and s > 0, then we have Vf. C HI, and for P G V^, 
we have the estimates 

\Piz)\ < ||P|| ( max \\z\\lA 

and for all z £ Ll, 



\\Xp{z)\\^,<2k{k-ir\\P\\ \\z\\^, ( max 



Moreover, for z and y in L], we have 

\\Xp{z)-Xp{y)\\^,<4k{k-ir\\P\\ ( max (\\y\\l, ,\\z\0)\\z - vW ■ 

n=l,...,k—2 ^3 

Eventually, for P G and Q G V^, then {P, Q} G P^_,_^_2 and 

\\{P,Q}\\ <2ke\\P\\ WQW 
and if P £ SV^ and Q G SV1, then {P, Q] G 5P^+^_2- 

Proof. The first estimate is a clear consequence of the definition of ||P|| . 
Assume that P is a homogeneous polynomial of degree k with coefficients air), 
T G Tfc. Assume that r G 2fc is fixed. The derivative of a given monomial 
z{t) = z{ti) ■ ■ ■ z{Tk) with respect to z{t) vanishes except if r C r. Assume for 
instance that t = = (d^)- Then the momentum zero condition implies that 
M{ti, . . . ,Tk-i) = — eC- Hence we can write with the previous notations 

dP 

(1 + Mr 



<^ll^ll E / (l + |r|)^|z(^,5)|de 



dz{T) 

Now as in (2.14) we have using the zero momentum condition 

(i + ir|)^<(i + i6i + ... + iefc-iir <(i-ir max (i + iCnir 

n=l,...,«;— 1 
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As in (2.15), we obtain the first estimate after integrating in r = (Ce)- The 
second estimate as well as the estimate on the Poisson bracket is proven similarly. 

■ 

We then consider operators Hq with continuous frequencies w(^). In the case 
of NLS we have uj{^) = l^p. 

Starting with this formalism, the formal analysis of section 3 is essentially the 
same with = othihUo) for some real function ah satisfying ah{x) ~ x for small 
X. Thus 

A^m = anihumm and Ao^(e) = -ah{hoj{0)i'{i)- 
and we set A(^) = ah{hu){^)). We consider the splitting method 

and we look for a real Hamiltonian function Z{t, (j), ■0) such that for all t < h we 
have 

¥p O = 

and Z(0) = Zq := ^o- 

As in section 3, plugging an Ansatz expansion Z{t) = J2e>i'^^ S^t the 
same recursive equations (3.6) to solve. The key lemma is the following: 

Lemma 6.4 Assume that 

Q{z) = [ a{T)z{T)dIk{T) 

is a homogeneous polynomial of order k. Then 

adzo(Q) = / iA{T)a{T)z{T)dIk{r) 

where for a multi-variable r = (ri, . . . , r,-) such that for i = 1, . . . ,r, Ti = {^i, 5i) S 
R'^ X {±1}, we have 

A(T) = 5lA(6) + --- + 5rA(Cr). 

Using this Lemma, we obtain that the condition (4.1) takes the form 

VA; = 3,...,r, ^relk, |A(r)| < 2^ - 5 

and the rest is extremely similar. For instance for the cubic NLS equation (6.1) 
with ah{x) = \//iarctan(x/\//i) and g{u,u) = \u\^ we get 
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Theorem 6.5 Let = {(j)o,(j)o) G be the initial datum in Fourier space 
(4)Q = uq) and define the sequence = {4>n,4>n) (the numerical solutions in 
Fourier space) by 



z 



n+l 



^'},o^\^{zn, n>0, 



where Aq is associated with the filter function ah{x) = arctan(2;/\//i) and P 
is associated with the nonlinearity j \u\'^ ■ Assume that for all n, = {4>n,4'n) 
remains bounded in . Then there exist constants ho and c such that for all 
h < ho, there exists a polynomial Hamiltonian Hh such that 

Hh{(l^nAn) = Hh{(t>A) + 0{eM-ch^^'^)) 

for nh < exp(c/i~^/^). Moreover, with the Hamiltonian H^^ given by 

H^h\^^^)= I ^arctan(|e|V^^)l</'(Ol'de+ 

i(|eiP + l6p-|r/ip-|r?2r 



X 



<A(6 ) <A(6 ) "0 ) (^2 ) dCi dC2d7?i d?72 



we have 



H^fl\(kn,4>n)=H^^\<t>,4>) + 0{h) 

over exponentially long time n < exp(c/i^^/^). 

Moreover, suppose that Uo G then there exists a constant C such that 

I \c\'\MO\'dC + ^ [ \MOM<c 

J|€|<Qh-i/4 ^/h J|5|>ah-i/'l 

over exponentially long time n < exp(c/i^^/^). 

Finally, assume that for all n, z" remains in B^j- with s > 3, then we have 

H{z"') = H{z^) + 0{h) 
over exponentially long time n < exp(c/i~^/^). 
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